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Abstract: We study in detail (p)reheating after multi-field inflation models with a partic- 
ular focus on N-flation. We consider a variety of different couplings between the inflatons 
and the matter sector, including both quartic and trilinear interactions with a light scalar 
^ I field. We show that the presence of multiple oscillating inflatons makes parametric reso- 

Q>^ I nance inefficient in the case of the quartic interactions. Moreover, perturbative processes 

do not permit a complete decay of the inflaton for this coupling. In order to recover the 
^ hot big bang, we must instead consider trihnear couphngs. In this case we show that 

^^ ■ strong nonperturbative preheating is possible via multi-field tachyonic resonance. In ad- 

^^ . dition, late-time perturbative effects do permit a complete decay of the condensate. We 

also study the production of gauge fields for several prototype couplings, finding similar 
results to the trilinear scalar coupling. During the course of our analysis we develop the 
rS I mathematical theory of the quasi-periodic Mathieu equation, the multi-field generalization 

C^ I of the Floquet theory familiar from preheating after single field infiation. We also elaborate 

on the theory of perturbative decays of a classical inflaton condensate, which is applicable 
in single-field models also. 
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1. Introduction 



Reheating after inflation is a crucial requirement for any model. Depending on how infla- 
tion ends and how the inflaton interacts, the transfer of energy from the inflaton sector into 
standard model (SM) degrees of freedom may proceed differently. The detailed physics of 
(p)reheating is rather model dependent and previous studies have uncovered a wide array 
of different quantum field theory (QFT) phenomena, including perturbative [l| and non- 
perturbative ||2|, |3|, ^ production of both bosons and fermions |p. At the nonperturbative 
level a rich variety of preheating processes are possible, including parametric resonance 
||2|, ^ ^, tachyonic preheating [^] and tachyonic resonance Q. Examples have also been 
studied where the dynamics of reheating involves distinctively string theoretic processes 
||9|, |10|, |ll[ and cannot be captured by ordinary QFT. (See also [12| for further studies of 
string theory reheating and [O] for other aspects of the endpoint of string theory inflation.) 

The requirement of successful reheating places a number of nontrivial constraints on 
inflationary model building. At the most basic level, one must ensure that the inflaton 
field completely decays into radiation after inflation, in order to match smoothly onto 
the usual hot big bang evolution. Moreover, one must verify that exotic relics (which 
might interfere with the successful predictions of big bang nucleosynthesis) are not over- 
produced. It is therefore crucial to understand the decay channels of the inflaton in order 
to assess the viability of any inflation model. Depending on the model, the process of 
preheating may also lead to various cosmological observables such as gravitational waves 
1 14] or nongaussianities |l^, p^ . 

Previous studies of (p)reheating after inflation have focused largely on the decay of a 
small number of fields in the inflaton sector Q. However, there are strong motivations to 
consider inflation models in which a huge number of fields play an important dynamical 
role. One might expect such models to arise naturally in string theory, where realistic 
compactifications often contain exponentially large numbers of moduli and axion fields 
1 17, 18, H^. Indeed, string theory appears to contain a rich landscape [20| and, in general, 
stringy inflation might involve the motion of a huge number of fields on a very complicated 
scalar potential |21, 22, ^, 24, 25]. Only in very special corners of this cosmic landscape 
would one expect inflation to take place along a single direction in field space. 

Multi-field inflation is also desirable from a phenomenological perspective. The "as- 
sisted inflation" mechanism (first discovered in ^^ and generalized to a much broader 
class of models in ]£^]) allows the collective motion of a large number of scalar fields (pi 
(i = 1, • • • , A^) to drive inflation, even when each individual field could not support inflation 
on its own. The essence of this mechanism is very simple. Suppose the fields are decoupled 
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Then each inflaton feels the Hubble friction induced by the collective potential energies 
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(1.2) 



but only the restoring force coming from its own potential 

'ci>i + 3H<Pi = -V^{^^) . (1.3) 

With sufficiently large A^ the motion of each (pi will be over-damped, even if each Vi{(j)i) is 
relatively steep, alleviating (to some extent) the fine-tuning problems associated with ob- 
taining flat scalar potentials (see, for example, pS]). This mechanism permits a realization 
of chaotic inflation in a regime where each (pi undergoes a sub-Planckian displacement in 
field space, whereas the collective radial excitation p (defined by p'^ = X^j'/'f) traverses a 
super-Planckian distance. Such a model goes a long way to ease concerns about the validity 



of the effective field theory description of large field inflation [29]. 



In |18] a particularly appealing realization of assisted inflation was proposed: N-flation. 
In this model one considers the dynamics of N axion fields (pi, each associated with a 
different broken shift symmetry (pi — )• (pi + Ci. In the region relevant for inflation, the 



potential takes the simple form (1.1) with Vi{(pi) — -^(pf- The potential of N-flation has 



been argued to be radiatively stable [|T^. (See also |31, 32, 33] for more details on the 
dynamics/phenomenology of this model.) N-flation is quite natural from the string theory 
perspective. In oriented critical string theories the massless 2- form field S^j^ can wrap a 
huge number of independent 2-cycles when the theory is compactified to four dimensions. 
Each such cycle leads to an axion at low energies. Explicit examples are known where the 
number of low energy axions may be as large as N = 0(10^) J^]. Using random matrix 
theory, it has been shown that the predictions of N-flation are largely independent of the 



details of the compactification ]19]. See, however, ]30] for complications associated with 
consistently embedding this framework into stabilized string theory vacua. 

In this paper we study (p)reheating after multi-field inflation. (See also J|^, ^].) 
We focus on N-flation as a specific example, however, many of our results will apply to 
multi-field models more generally. We proceed phenomenologically, considering a variety 
of different couplings between the inflatons and matter. For simplicity, we first consider 
reheating into a scalar field x through prototype interactions of the form 

TV 2 ^ 

These two kinds of couplings yield very different results for the decay of the inflaton, both 
in the perturbative and nonperturbative regimes. In the case of 4-leg interactions, nonper- 
turbative preheating processes are suppressed by the de-phasing of the N-flaton oscillations. 
We show that reheating must therefore proceed via perturbative decays for this choice of 
coupling. This is in contrast to the single-field case, were nonperturbative parametric res- 
onance may occur. On the other hand, we show that the 3-leg interactions do allow for 
violent nonperturbative particle production via tachyonic resonance. The presence of a 
huge number of oscillating fields in the inflationary sector does not prohibit preheating via 
such couplings. During the course of our analysis we develop the mathematical theory of the 
quasi-periodic Mathieu equation, the generalization of the usual Floquet theory to multi- 
field inflationary models. We study quantitatively the structure of the stability/instability 



charts for multi-field inflation and comment on the dissolution of the stability bands in the 
presence of multiple incommensurate frequencies. 

Having ruled out nonperturbative preheating effects in the case of the 4- leg interactions, 
we next study the perturbative QFT decays of the coherent oscillating N-flatons. We find 
that the 4-leg interaction does not permit a complete decay; some relic density of inflaton 
particles always freezes out and the late-time universe is typically cold and devoid of life 
as we know it. This result is completely analogous to the single-field case, which we re- 
consider. 

Next, we study the perturbative decays of the coherent inflaton oscillations via 3-leg 
interactions, which will be relevant during the final stages of reheating and also when the 
inflaton oscillations have small amplitude. Here we show that the inflaton condensate is 
permitted to decay completely and the late-time universe may be compatible with the usual 
hot big bang scenario. Again, this result is completely analogous to the single-field case. 

Finally, we also consider coupling the N-flaton particles to gauge fields. We consider 
the following prototype interactions between the infiatons (pi and a U(l) gauge field A^^ 
with field strength F^j^ = (9^j4jy — d^A^: 

N , N 

Here F^j^^ = ^^^'^'"^Fpa is the dual field strength tensor, e^'^'"^ is the anti-symmetric Levi- 
Civita symbol and the Mj's are symmetry breaking scales. The former axion-type coupling 
to F A -F is directly relevant for N-fiation. In this case we find that the production of 
gauge fields is very similar to the 3-legs scalar interaction above; the inflaton is permitted 
to decay completely and nonperturbative preheating is possible via tachyonic resonance, 
depending on the model parameters Mj. For completeness, we also consider the coupling 
^j (j)iF'^ /Mi, which might be relevant for string theory models in which inflation is realized 
via the motion of a large number of moduli flelds. 

The organization of this paper is as follows. In section we review the N-fiation 
model and its string theoretic realizations. In section we discuss the possibility of non- 
perturbative preheating after N-fiation. During the course of our analysis, we develop the 
mathematical theory of the quasi-periodic Mathieu equation, a multi-field generalization 
of the Floquet analysis that is familiar in the preheating literature. Next, we turn our 
attention to perturbative decays. In section Q we review the perturbative decays of the 
coherent oscillations of the infiaton in single field inflation and compare this to the more 
familiar case of decaying infiaton particles. In section ^ we discuss the perturbative de- 
cay of the N-flatons via 4-leg interactions, showing that the infiaton is not permitted to 
decay completely. In section |6| we consider instead perturbative 3-leg interactions show- 
ing that, this time, a complete decay is possible. In section |^ we study the production 
of standard model gauge fields via both perturbative and nonperturbative processes after 
N-flation. In section |8| we discuss some possible applications of our results outside of the 
context of N-fiation and also study the production of gauge fields via "modular" couplings 
to the infiatons. Finally, in section ^ we conclude and speculate on directions for future 
investigation. 



2. The N-flation Model 

2.1 Homogeneous Dynamics 

In this section we briefly review the N-flation model and its stringy reahzation. The idea 
of p^] was to reahze multi-fleld assisted inflaton using axions. To flrst approximation, the 
potential takes the separable form ( |l.lD where the Vi are periodic and arise solely from 
nonperturbative effects: 



VM,) = A: 
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Here fi is the axion decay constant and Aj is some dynamically generated scale. For Aj 
well below the ultra-violet (UV) scale we can neglect higher harmonics in the potential and 
also cross terms coming from multi-instanton effects, which might spoil the separable form 
( |1 . 1|) . For small fleld values (pi <C fi we can expand the potential Vi as 



m: 
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V^i<P^) = ^cPj, (2.2) 

where m,j = 27rAf/fi. This potential provides a realization of chaotic inflation for the 



collective variable p = \ j^Zi ^li which may experience a super-Plan ckian displacement in 
field space, even when the individual field displacements are sub-Planckian. 

The homogeneous dynamics of the N-flaton fields are described by the usual Klein- 
Gordon equation 

4>i{t) + m4>i{t) + mj(t)i{t) = 0. (2.3) 

During the slow roll phase the friction term dominates and one can easily see [^] 

Mt) _(<hiJ^Y^"^\ (2.4) 



4>i{tQ) \<t>j{to 

This shows that even if the fields begin near the same point in field space they quickly 
spread out during infiation: the heaviest fields damp rapidly to the minimum and drop 
out of the inflationary dynamics while the lightest fields remain frozen until the end of 
inflation. We illustrate these dynamics in Fig. |l|. 

Once inflation has ended, the lightest N-flaton fields (those which have not been 
damped to zero during infiation) will oscillate about the global minimum of the poten- 

2 

tial V = ^j ~2^4'i- It is this phase of post-infiationary oscillations which is relevant for 
(p)reheating. To describe these oscillations, note that equation ( |2.3D can be written as 
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Since m? 3> H^ , H during preheating the surviving N-fiatons will oscillate as 
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Figure 1: Evolution of several N-flaton fields showing both the inflationary period and post-inflation 
oscillations. This figure illustrates the damping of heavier fields during inflation. At the onset of the 
oscillatory phase relevant for (p)reheating, only the lightest fields remain dynamical. Notice that 
these are nearly degenerate in mass, and therefore their post-inflationary oscillations are nearly in 
phase. For illustration we take N = 1500 fields with equal energy initial conditions at the beginning 
of inflation. 



similarly to the single field case |^]. The coefficients (/)o,t and phases 9i are determined 
by the initial conditions at the end of N-flation. The scale factor, averaged over several 
oscillations, grows as o(t) ~ t"^'^ while the energy density in the N-flaton fields decreases 



as Eii 



E^ 
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~ a~ . Hence, the coherent oscillations of the inflatons after 
N-flation behave very much as massive, non-relativistic particles. 

The homogeneous dynamics of the N-flatons, including the phase of post-inflationary 



oscillations, were studied in detail in |34|. For reasonable choices of mass spectrum and 
initial conditions (at the beginning of inflation) it was shown that, at the onset of the phase 
of oscillations, more than 90 percent of the energy is confined to the lightest 10 percent 
of the fields in a very narrow mass range. The remaining heavier fields will have already 
settled down at the minimum by this stage. Since the fields which survive to the end of 
inflation are nearly degenerate in mass, their post-inflationary oscillations are nearly in 
phase; see Fig. |l|. It is this subset of fields that plays an important dynamical role during 
(p)reheating. More on this later. 

2.2 String Theory Realization 



In |18] it was argued that the potential (2.2) is radiatively stable. However, since the 
flatness of the inflaton potential may be sensitive to Planck-suppressed corrections, it is 
useful to comment on how this scenario might arise in string theory. (See also |^0| and [^].) 
For the sake of being concrete, suppose we consider type IIB string theory and stabilize 
the moduli using the method of KKLT |39|. (Similar results are expected in the stable 
HA vacua of [40|.) In the dimensional reduction, there arise A^ = hi^i axions (pi associated 
with the 4-cycles of the Calabi-Yau manifold.^ In a supersymmetric compactification, 
these axions combine with the volumes of the 4-cycles tj to form the complex parameter 
Ti = Ti — icpi. The KKLT super-potential isW = W^ + ^^^i ^jC""'-^' where Wq, Ai depend on 



""^The meaning of the tilde notation cpi will become clear shortly. 



the dilaton and complex structure. The exponential terms in the super-potential arise due 
to nonperturbative effects (one for each 4-cycle). To leading order, the Kahler potential 
can be written as K = — 21n(V) where the Calabi-Yau volume V depends on the 4-cycle 
moduli Ti in a non-trivial way. For small values of the axions (pi, the action takes the 
schematic form 
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(2.7) 



In general, neither the mass matrix {M'^)ij nor the metric on field space Gij is diagonal. 



To put (2.7) into canonical form we must perform a rotation in field space. In |19| this 
diagonalization was performed explicitly. For our purposes, it suffices simply to note that 
one must generically perform a transformation of the form 

4>i = Y^af(t>j (2.8) 



to put ( |2.7| ) in the canonical form with N-fiaton potentials given by ( |2.2| ). In |19] it 
was shown that the masses rrii are distributed according to the Marcenko-Pastur law, 
independently of the microscopic details of the compactification. This distribution, which 
we assume throughout this work, depends on a single parameter (3 which counts the number 
of axions divided by the dimension of the moduli space. Where explicit numerical values 
are required, we assume the value f3 = 1/2, which was advocated in |19| |. 

An important observation is that the interactions between the N-flatons and visible 
matter will typically be diagonal in the basis (pi which is associated with the geometrical 
4-cycles of the Calabi-Yau manifold. For example, we expect a coupling of the form (pjFAF 
to gauge fields living on a D7-brane wrapping the j'-th cycle of the Calabi-Yau manifold. 
As a result of the field rotation ( |2.8| ), we generically expect that all canonical N-flatons 
couple to each visible matter sector. 



3. Preheating After N-flation 

3.1 Couplings to Scalar Field Matter 

In this section we study the possibility of strong nonperturbative preheating effects after N- 
fiation. To proceed, we must first couple the inflatons to some matter fields. For simplicity, 
we first model matter by a single light scalar field x- (Later on we will consider more realistic 
couplings to brane-bound gauge fields.) Since we do not expect the interaction basis (pi 
and the mass basis (pi to be aligned, we should consider couplings between x and all the 
N-flatons. We therefore propose a prototype Lagrangian of the form 
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(3.1) 



The term Ax'^ is required to keep the potential bounded below, but otherwise will play 
no role in our discussion. Treating the N-flatons as classical background fields ( |2.6| ), the 
equation of motion for the linear quantum fluctuations of the x field takes the form 



V2 

X + ^Hx + -^X ■ 



Y.{g}^\{t)-rCJi^^{t)\x = ^■ 



(3.2) 



We introduce the variable X(t,x) = a^'^(t)x(i,x) and write the equation of motion (| 
in the form 



Xfc(i) + 
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where A- •" i'^" ' u2\ _ 3 P 

Since the N-flaton oscillations ( |2.6| ) behave as non-relativistic matter, P 

ignore the quantity A in (|3.3|). The effective mass for the Xk modes is 



2| + H^^ = j^7 with P the pressure of the dominant inflaton background. 

and we can 
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(3.4) 



If the N-flaton masses were precisely degenerate, rrii = m, then the oscillatory contributions 
to (3^) would add coherently, mimicking the well-studied single- field case |3|, g]. More 



realistically, however, the inflaton masses m-j should be sampled from some distribution such 
as the Marcenko-Pastur law [19|. In this case the time-dependent behaviour of the effective 



mass (3.4) may be quite different from the single-field case and, therefore, preheating may 
proceed differently. To understand the behaviour of the solutions of ( |3.3| ) we should first 
study Mgg(i) for parameter choices that are reasonable for N-flation. In passing, notice 



that the summation over N-flaton field indices in ( |3.3D and (3.4) need only include the 
lightest N-flatons which remain dynamical until the onset of post-inflationary oscillations. 



In constructing solutions of (3^) there are two limiting cases which may be of interest, 
depending on which interaction in (|3.l| ) dominates. It may happen that the 4-leg interac- 



tions 'Ylid'i'i^'iX^ dominate over the 3-leg interactions J^i^i'l^iX^^ ^^ '^ise versa. Let us first 
study the behaviour of M^g{t) in each limiting case, to demonstrate the absence of strong 
parametric resonance preheating via the 4-leg interactions and the possibility of strong 
tachyonic resonance preheating via the 3-leg interactions. Later on, we will confirm these 
results by developing the mathematical theory of the quasi-periodic Mathieu equation. In 
that case, we will see that both limiting cases can be treated simultaneously by means of a 
simple device. 

3.2 Four-Leg Interactions: Absence of Strong Adiabaticity Violations 

Let us first study the time dependence of the effective mass (^^) in the limit where the 
3-leg interactions are vanishing, that is cjj = 0. Now (|3.4| ) takes the form 



2^2 
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M^,{t) = 5: gf^it) = E 5^ «i-' (-^* + ^^ 



(3.5) 



In Fig. g we illustrate numerically the typical time-dependence of the effective mass (^ 
for A'^ = 600 inflatons with the Marcenko-Pastur mass distribution [19| and randomly 
distributed initial conditions (/>i(to) at the beginning of N-flation (the qualitative behaviour 
is independent of this specific choice). 




Figure 2: Evolution of the effective mass 
^4ff(0 - Eaf^t^Kt) of the X field for N ^ 
600, assuming randomly distributed initial 
field values at the start of inflation. Only the 
oscillatory part of the motion relevant for re- 
heating has been illustrated. For simplicity, 
we have set all the coupling constants g^ = 1. 
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Figure 3: The time evolution of the quantity 
Mes/M^ff for the effective mass M^g{t) = 
'^i9i4'i{t) ^-nd the same parameters as 
Fig. 1. This quantity is a measure of adi- 
abaticity of the mode functions. For param- 
eter choices relevant for N-flation strong vi- 
olations of adiabaticity are absent. 



The spread of the masses rrii in a realistic N-flation model prevents the effective mass 
( |3.5| ) from crossing zero and impedes the violations of adiabaticity that are associated with 
broad band parametric resonance. To illustrate this we plot the quantity Mes/M^g, which 
measures the adiabaticity of the mode functions, in Fig. ^. We have verified the smallness 
of the adiabaticity parameter in the large- A'^ limit for several choices of initial conditions 
(see also |^^). Therefore, we do not expect any strong nonperturbative preheating effects 
for the 4-leg interactions 'Ylidf'i^'iX^^ assuming A'^ ^ 1 and the Marcenko-Pastur mass 
spectrum. This expectation was verified numerically in [^^. Below, we confirm these 
results using the mathematical theory of the quasi-periodic Mathieu equation. 

3.3 Three-Leg Interactions: Possibility of Tachyonic Resonance 

Let us now consider the opposite regime, wherein the 3-leg interactions in ( p.lD dominate. 
We therefore seek to solve the mode equation (|3.3| ) with gf = 0. In this limiting case the 
effective mass (|3.4|) takes the form 



^eff(*) = E ^*<^^(^) - E ^ ^i^ ("^^^ + ^^ 



The resulting mode equation is 
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Xk{t) = 0., 



(3.6) 



(3.7) 



where A = |tit — can be neglected. The behaviour of the effective mass ( |3.6D is 
very different from the 4-legs case, equation ( |3.5D . In the presence of dominant 3-leg 



interactions the effective mass-squared can go tachyonic, Mgg(i) < 0, at certain stages 
during the N-flaton oscillations. In the single field case, this leads to the well-studied 
tachyonic resonance preheating Q. During a tachyonic phase the fluctuations Xk{t) will 
display exponential growth, which is associated with the explosive production of x particles 
by spinodal instability. 

The strength of the nonperturbative tachyonic resonance effects associated with ( |3.7] ) 
depend sensitively on the couplings ai and the initial conditions at the onset of reheating. 
As discussed previously, the heaviest N-flatons drop out of the dynamics rapidly during 
N-flation. During the post-inflationary oscillations, the majority of the energy is carried 
by the lightest 10 percent of the N-flatons with a very narrow range of masses. Since these 
N-flatons have nearly degenerate mass, they start to oscillate at roughly the same time, 
leading to a picture where the oscillations ( |2.6| ) are more-or-less in phase, as illustrated 
in Fig. ||. This approximate coherence results in a time-dependence for the effective mass 
( |3.6| ) which is nearly identical to the results one would obtain for a single field inflation 
model. ^ The resulting dynamics of the effective mass (^^) are illustrated in Figure 0. 




Figure 4: The post-inflationary time-dependence of the effective mass M^g(t) = J2i o'i4>i(.'t) ^or the 
3-leg interaction. For simplicity, we assume ai = a. This figure shows that the post-inflationary 
oscillations in M^g(i) are qualitatively identical to what one would obtain for single- field inflation 
(corresponding to the sin(TOf) solution), aside from a slowly time-dependent frequency and more 
rapid damping of the amplitude. In plotting the analytic solution, we have numerically matched 
tlie amplitudes at the first mimima and computed fh from the first half oscillation of the multi-field 
solution. 



From Fig. ^ we see that the effective mass ( p.6| ) looks like that of a single field with 
a slowly time- varying frequency and an amplitude damping faster than would be expected 
from the expansion alone. Both of these effects are due to the de-phasing of the various 
N-fiaton fields (this may be verified analytically also). 

In summary: we have seen that a realistic choice of N-flationary parameters leads to 
a behaviour of the effective mass (|3.6D in the case of 3-legs interactions which is nearly 



"As we saw above, in the case of the 4-leg interactions (3.5) any slight de-phasing of the N-flaton 



oscillations makes violations of adiabaticity exceedingly unlikely in the limit A'' ^ 1. This occurs because 



(3.5) is a sum of N positive definite terms, when N is large zero-crossings become very rare. However, in the 



case of the 3-leg interactions the effective mass (S.C) is very different; it is a sum of A'' terms with indefinite 
sign. Hence, small de-phasing of the N-flaton oscillations do not conspire to forbid tachyonic regions. 
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identical to the usual single- field result. Hence, preheating of the X^ fluctuations will 
proceed via the same tachyonic resonance mechanism which was studied in detail in p, 57 1. 



We do not repeat this analysis here. On the other hand, we can imagine a scenario wherein 
a coupling of the form ( |3.6D is relevant for preheating, however, the various inflaton fields 
have amplitudes which are more-or-less randomly distributed at the onset of reheating.^ 
For example, such a situation might arise in multi-field models where inflation proceeds via 
some extremely complicated trajectory through the string landscape |21, |2^, ^, ^, pSJ . 



Therefore, let us briefly consider a toy scenario (not necessarily relevant for N-flation) where 
N = 5 inflatons have random initial displacements satisfying 0.2Mp < |(/>o,j| < Mp at the 
onset of reheating. We also assume a random mass spectrum. Our choice of parameters for 
this example is given explicitly in Table |l[ By choosing a somewhat large hierarchy between 
the heaviest and lightest field we can obtain results which are qualitatively different from 
the single-field case. 

Table 1: Masses and Initial Displacements of Inflatons 



Mp 



1 2.43 0.533 

2 5.63 0.412 

3 9.88 0.739 

4 15.85 0.544 

5 26.67 0.701 



In Fig. ^ we display the post-inflationary oscillations of the heaviest and lightest infla- 
ton fields for the example parameters in Table |l[ In Fig. ^ we plot the time-dependence of 
the effective mass-squared ( |3.6| ) for this same example, showing the presence of tachyonic 
regions where Mgg(t). These regions leads to violent particle production via tachyonic 
resonance. From Fig. M it is also clear the behaviour of M^g{t) is qualitatively different 
from the single field case (compare with Fig. ^). 

The tachyonic regions illustrated in Fig. g and the non-periodic oscillations of the ef- 
fective frequency lead to unstable growth of the fluctuations Xi^(t) by tachyonic resonance. 
The exponential instability of the Xk inhomogeneities is associated with explosive particle 
production. In Fig. and Fig. |^ we have illustrated the unstable behaviour of the solutions 
Xfc(t) of ( ^ ) taking, again, the parameter values in Table ||. These figures illustrate the 
rapid growth of inhomogeneities, showing that this new type of multi- field tachyonic reso- 
nance is very efficient. Within a very short time, back-reaction effects become important 
and one must turn to lattice simulations to fully capture the dynamics. 

3.4 The Quasi-Periodic Mathieu Equation 

We have argued that strong nonperturbative effects are absent for the 4-leg interactions 
( |3.5| ) , whereas tachyonic resonance may proceed in the case of the 3- leg interactions (^^) . 



^This is, of course, very different from assuming randomly distributed inflaton field values at tlie onset 
of inflation. 
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Figure 5: The evolution of the most mas- 
sive (solid) and least massive (dashed) back- 
ground inflatons for five fields for the param- 
eters given in Table 1. As we see, each infla- 
ton behaves as a damped harmonic oscillator 
as expected. As usual, the fields values are 
in units of Mp and the time is z = rriaunt- 




Figure 6: The evolution of the effective mass 
M^ff = ^^ (Ji(j)i{t) of the X field for the ex- 
ample in Table 1. We assume that all cou- 
plings are equal: Ui = a for all i. This fig- 
ure illustrate the presence of multiple tachy- 
onic regions during the inflaton oscillations. 
The dynamics of the effective mass is very far 
from periodic, which implies that this exam- 
ple is qualitatively different from the single- 
field case. 




Figure 7: The unstable growth of the quan- 
tum field Xk(i) due to multi- field tachyonic 
resonance effects. For illustration we take 

k = 40r71min. 




Figure 8: The unstable growth of the quan- 
tum field X\s_{t) due to multi-field tachyonic 
resonance effects. For illustration we take 

k = 9.608TOmin. 



We will now re-consider these results using the mathematical theory of the quasi-periodic 
Mathieu equation ||4^, 5C, 47, ^J. Using this perspective, we can unify the discussion of 
both the 4-legs and 3-legs interactions in ( |3.lD . To see this, note that, in the limit where 



the 4-leg interactions dominate, equation (3.3) yields 
Xk{t) + 



— + 2_^ Ts77^ ^^"^ V"^j* + 0i) + i\ 

i 



a\t) 



^fc(t)=0. 



(3.8) 



On the other hand, in the limit where the 3-legs interactions dominate, we would have 



Xk{t) + 



— + 2_^ :k/'),,\^ sm [rriit + Oi) + I\ 



a 



a 



^l^(t) 



Xk{t) = Q. 



(3.9) 



In both cases, if we ignore the expansion of the universe, an appropriate redefinition of our 
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time parameter puts the equation in the form 



,^,X.(.) + 






j4fc — 2 > qi cos 2 — z + 

Z — J \ rn 



m 

i=l 



Xfc(z) = 0, (3.10) 



where m is some convenient mass scale (such as the mass of the heaviest or hghtest inflaton, 
or the average mass of the infiatons) and A'^osc is the number of N-flatons which play 
a significant role in reheating (that is, A'osc counts only those lightest N-flatons which 
oscillate at the end of inflation). Therefore, we can present of unified treatment of both 
limiting cases by developing the theory of equation ( p. 101 ) . 

Let us now discuss the solutions of the quasi-periodic Mathieu equation ( |3.10| ). In 
the simplest case of a single infiaton field with rrii = m and qi = q. Then, neglecting the 



expansion of the universe, equation ( 3.10 ) reduces to the Mathieu equation [^] which is 
well-known from studies of preheating after inflation y] . This equation is known to posess 
instability bands, for certain values of A^. and q the solutions display exponential growth 

Xk{z) = €>'>'' Pk{z), (3.11) 

where Pk{z) is some periodic function and the real number //^ is the Floquet exponent 
which determines the rate of the instability. The stability chart of the Mathieu equation, 
depicted in the left panel of Figure |9|, can be used to explain the resonance structure of 
preheating for single fleld inflation [g, y] . The shaded regions of this chart display parameter 



values for which the solutions of ( 3.10 ) have Floquet exponent /i^ 7^ and hence display 
unstable growth. 

For the case at hand, multiple inflaton flelds, there are two cases to distinguish de- 
pending on the values of the oscillation frequencies rrii/m: either these are all rationally 
related or else at least two of them are incommensurate. In the flrst case, the effective 
frequency for Xk is periodic and Floquet's theorem ensures that for certain values of the 
parameters Aj. and qi the solutions will again grow as ( p.ll| ). Instead, if two of the masses 
are irrationally related, then the effective frequency is no longer periodic and Floquet's 
theorem no longer applies. When this occurs, ( |3.10|) is known as the quasi-periodic (QP) 
Mathieu equation. We will study the theory of this equation below. With some abuse of 



language, we will refer to ( 3.10| ) as the QP Mathieu equation and the exponential growth 



rates as the Floquet exponents regardless of the actual mass spectrum. One may worry 
that the growth is no longer exponential in the QP case. However, since an infinitesimal 
change in the masses can change their ratios from commensurate to incommensurate and 
vice versa, we expect exponential growth in this case also. 

In the QP case (incommensurate frequencies), we can also make a further distinction: 
all of the masses and their sums and differences are incommensurate, or at least one pair of 
masses (or some pair of sums or differences) are rationally related. The second of these can 
lead to the phenomenon of resonance trapping |^^ and may play a role when both 3-legs 
and 4-legs interactions are relevant during preheating. A detailed study of these effects is 
beyond the scope of this paper and we will instead concentrate on cases where only a single 
type of interaction is relevant during preheating. 
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So far, we have ignored the expansion of the universe. Adding the expansion makes the 
resonance parameters Ak and g, time dependent (see below), and the results of Floquet's 
theorem will no longer apply, even for rationally related masses. However, notice that 
H <^ rui during reheating, thus the time scale for variation of a{t) is much longer than 
the period of oscillations in the effective frequency. We can therefore employ an adiabatic 
approximation wherein we solve ( 3.10| ) using the standard techniques but allow A/., qi to 



vary slowly in the solution as a result of the expansion of the universe. This leads to a 
picture where the dynamics of a{t) causes the resonance parameters qi to "flow" across the 



stability /instability chart. For a similar discussion in the single field case see [56|. More 
details on this shortly. 

3.5 Stability Chart and Floquet Exponents 

Having given the general picture of the relation between the QP Mathieu equation and 
reheating, we now turn to develop a mathematical understanding of the stability/instability 
chart associated with ( 3.1[1| ) in the interesting multi-field case. Later on, we will apply these 



results to the specific reheating models under consideration. 

In the case A'osc > 1 the equation ( |3.10| ) is immensely more complicated than the 
usual Mathieu equation. As a warm-up, we first consider the illustrative case A'osc = 2 
with qi = q2 = q and neglect the phases 6i =62 = 0. (This very restrictive case may 
actually have some relevance for N-flation since only a relatively small number of N-flatons 
with nearly degenerate masses are expected to participate in preheating.) For this simple 
example we have solved ( |3.1[1| ) numerically for a wide range of parameters Ak and g, taking 
mi/m = 1 and m^jm, = 0.3475. We have numerically constructed the stability/instability 
chart, and the result is displayed in the right panel of Figure |9|. 

As can be seen from Fig. ^ the non-commensurate "extra" oscillatory contribution to 



(|3.10D leads to the formation of "new" instability regions at small g^j, called Arnold tongues. 



These emerge from the A-axis with locations given by 



A 



Errii 
rii — 
m 



(3.12) 



where each ni is an integer and the sum runs over the oscillating inflatons ^, ^]. We 
will refer to ^^ |nj| as the order of the resonance. For qi <^ 1, the tongue labeled by the 
vector n = (ni, ■■■,nN^^^) has width ~ Yli \qi\ ■ If the masses are all rationally related, 
then eventually the higher order resonances will be located at the same A value as a lower 
order resonance and the number of Arnold tongues will be finite. However, if we have 
incommensurate masses, then the tips of these instability tongues will densely fill the A 
axis. Thus, one might naively expect that in this case all of the modes will be unstable 
and the resonance will be very efficient. However, even if the stability bands completely 
dissolve, this does not necessarily mean that preheating will become significantly more 
efficient. It might be that the Floquet exponents fi in the dissolved regions are very small. 
Let us now discuss the Floquet exponents /i associated with the quasi-periodic Math- 
ieu equation. In order for the exponential growth ( f3.ll ) to be important, the Floquet 



14 





Figure 9: In the left panel we show the stability/instability chart for the usual Mathieu equation, 
that is ( 3.10 ) with TVosc = 1, m-i = m and qi = q. In the right panel, we display the analogous 
stability/instability chart for the quasi-periodic Mathieu equation, that is (5.1C) with A*"osc = 2, 
m2 = 0.3475mi and qi ^ q2 = q. In both panels, the white regions of parameter space correspond 
to stable oscillatory solutions, whereas the shaded regions correspond to parameter values for which 



the solutions of ( 3.10 ) display unstable, exponential growth. We see that the inclusion of an 
additional oscillatory contribution to the effective frequency of equation ( 3.1C ) leads to the formation 
of instability "tongues" which permeate the regions that were formerly associated with stability. 
There are, in fact, an infinite number of such instability tongues suggesting that the stability region 
might be dissolved. However, we only display a finite number of instability tongues, neglecting 
higher order effects which lead to very weak growth that cannot be accurately captured by finite- 
duration numerical simulations. The relevance of the lines A — 2q, A — Aq and A ^ Qqis discussed 
below. 



exponent /i must be large enough that the growth e'^^ is appreciable over the time of re- 
heating. Moreover, the expansion of the universe will cause modes to "flow" through the 
stability/instability chart, weakening parametric resonance effects as a{t) grows. If// is not 
sufficiently large, then there will be no important growth of the mode before it is red-shifted 
out of the instability band. Let us consider the behaviour of the Floquet exponent in our 
two-oscillator example with A'^osc = 2 with qi = q2 = 0.05 and m2 = 0.3475mi. Fig. 10 
shows ^ as a function of A. The two first order resonances (^^ \ni\ = 1) are much stronger 
than the second order (^^ |nj| = 2) resonances, while the third order resonances are too 
weak to be visible. Hence, in this case the fact that the Arnold tongues densely fill the A 
axis in Fig. |9| turns out to be irrelevant: the "instability" associated with the higher order 
tongues is far too weak to be important for reheating. 

It is interesting to compare the Floquet exponents for our example quasi-periodic 
oscillator to the usual Mathieu equation, in the regime q <C 1. From Figure |l^ we see that 
for the resonances which are common to both the QP Mathieu and the periodic Mathieu 
equations {i.e. those with only a single n^ non-zero), the Floquet exponents are nearly 
identical for both equations. The effect of the multiple frequencies is simply to create new 
resonance bands which do not overlap with those present in the periodic Mathieu equation. 
Since these resonances are all of order 2 or greater, they are much weaker than the first 
order resonances and are thus sub-dominant. 

Thus far, we have focused on the Floquet exponents in the regime g^ ^ 1, which is 
associated with weak resonance. In the single field case, g ^ 1 would be the condition for 
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Figure 10: Floquet exponents for the case of qi = q2 = 0.05 and a range of A values. The locations 
of the instability bands are clear, with the two widest bands corresponding to the two first order 
resonances. The remaining three resonances are all second order and correspond, from left to right, 
to the (711,^2) = (1,-1), (0,2) and (1,1) resonances, where the first index is for the mi — m 
oscillator and the second for the m2 = Q.3475m oscillator. Also shown are the Floquet exponents 
for the case of a single oscillator (corresponding to the periodic Mathieu equation) with a mass 
TOi = m (dashed green curve) and TO2 = 0.3475?7i (solid red curve). 



a strong parametric resonance. How should this condition be generalized to the multi-field 
case? If we denote the Floquet exponents for the periodic Mathieu equation with mi/m = u 
as fJ-uiiA, qi) then a simple re-definition 2: — )• wz shows that Hij{A, q) = Hi^=i{A/uj'^,q/u'^)iLi. 
This demonstrates that the effective resonance parameter (when determining Floquet ex- 
ponents) for an oscillator with mass m, is given by qf = qi {m/rrii) . Hence, we expect 
(roughly) that strong resonance will occur in instability tongues involving the i-th oscillator 
(those with rii 7^ 0) whenever qf ^ 1. However, as we explain below, this intuition may 
need to be modified when q becomes sufficiently large. 

Consider now the approach to the strong-resonance regime of our toy two-oscillator 
example. From the stability diagram, we see that increasing q causes the various resonance 
tongues to open up and begin to overlap. This is a qualitatively new feature relative to 
the periodic case (where the tongues do not overlap). We have illustrated the behaviour of 



the Floquet exponents in this overlap region for our two-oscillator example in Fig. 11. (We 
compare this two-oscillator case to the usual Mathieu equation with m = 1 01 m = 0.3475.) 
Prior to the overlap, the strength of the resonances can still be estimated by considering 
the parameter g|^ = qj(m/mj)^ relevant for the single oscillator case. However, once 
multiple oscillators contribute to the instability for a fixed value of A and q, the Floquet 
exponents become complicated functions of A (for a fixed q). When this occurs, it is much 
more difficult to determine the effective resonance parameter controlling the strength of 
the instability. 



Let us estimate the qi values at which the overlap displayed in Fig. 11 occurs. Clearly, 
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Figure 11: Floquet exponents for several different equations as a function of A for a range of q 
parameters. In the upper panel are the exponents for the QP Mathieu equation with qi — q2 = q 
and the same masses as in Figure pi The bottom left panel is for the periodic Mathieu equation 
with m — 0.3475. The bottom right panel is for the periodic Mathieu equation with m = 1. 



if we have q'[^ > 1 for one or more of the oscillators, then there will be many "open" 
resonance bands and we would expect them to overlap with each other. For our illustrative 
2 oscillator model, this is what occurs. For a large number of oscillators, the precise 
location of the overlap region is very difficult to estimate if we have g| > 1. Therefore, we 
will consider instead the limit when q*^^ <C 1 for each of the oscillators when the overlap 
occurs. Since the most important resonances are the first order ones in this case, we will 
consider where the first order resonances overlap. For qf <C 1, the boundaries of the ith. 
stability tongue are given by j4min = m^/m? — qi and A^ax = mf /m? + q^. [|^, 50] Define 



-^^l- 



m- — mf. Then the ith and {i + l)th first order bands overlap when 



+ ^ 



li+l 






(3.13) 



For the purposes of illustration, suppose that the mass splittings and qis are all equal 
Denote drnf^^i = 5m? ai 
bands will overlap when 



Denote Smf^^i = 5m? and qi = q. Then the above tells us that the first order resonance 



5m? 
2m? 
The effective q value in the ith first order band is then 



eff, overlap 



2^ 



(3.14) 



(3.15) 
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Figure 12: Floquet exponents for q — 200 for the QP Mathieu equation and the Mathieu equation 
with two different oscillation frequencies. In the QP case, we see that far fewer modes are stable 
and that the resonance remains strong for larger values of A. 



Hence, the above estimate is consistent provided the mass distribution is not too broad 
(i.e. 5m? <C mf for each oscillator). 

Finally, consider increasing q even further, into the regime where the resonance bands 
are no longer in the process of overlapping. Figure 12 demonstrates that the Floquet 



exponents become much simpler functions of A and the complicated effects of interference 
wash out. Relative to the Mathieu equation, we see that the resonance is more efficient 
for essentially all choices of A. Additionally, there are fewer stable modes and that the 
resonance extends over a larger range of A values. 

Before applying the above mathematical framework to the physical problems at hand, 
one final remark is in order. Our above analysis of the Floquet exponents took a global 
viewpoint in that they were computed by determining the growth of the mode functions 
over long time intervals. However, during reheating, the resonance parameters evolve 
and only the early parts of the evolution of the mode functions tend to be relevant. As 
well, when the resonance is strong, the mode functions will grow very quickly and our 
linear analysis soon ceases to be valid. Again, it is only the early time evolution of the 
mode functions that tends to be relevant. Thus, some knowledge of the form of the mode 
functions themselves (not just their average growth) is required in order to truly apply the 
above formalism to reheating. However, here we will be content primarily with deciding 
whether or not preheating is possible, in which case a knowledge of the asymptotic Floquet 
exponents should be sufficient. 

3.6 Thee-Legs Interaction: Geometrical Interpretation 

Having developed the mathematical theory of the QP Mathieu equation, let us now ap- 
ply this insight to study reheating after N-flation. We begin by considering the 3-legs 
interaction in the context discussed above. Equation ( |3.9[ ) takes the form ( p. 10 ) with 
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Figure 13: The paths through the instabihty chart traced out by several mode lunctions (with 
different values of k) for the case of a 3-legs couphng between the inflatons and scalar field. We 
have taken our initial value of q to be qo — 2. As time passes, all of the modes move towards the 
origin. 



2z = mt + 7r/2 and the resonance parameters (ignoring terms involving H/m^ which is the 
condition for A and qi to evolve adiabatically) 



9 9 ' 



lOi 



50,j 



m?a?i'^ 



(3.16) 



Recall our picture of the parameters of our mode functions tracing out curves in the 
(A, q) parameter space. In this case, the initial parameter values are distributed along 
the vertical line qi^ = '^^ ^jl with ^ > 0. Subsequently, each mode will evolve along 

the curve A^ = A, 



2 3/2 
4/3 



ifc,o I "^ ) • This evolution is illustrated in Fig. 13 for a particular 
choice of initial conditions qi^ (determined by (j)Q^i and (Ji) and several k values. As can 
be seen from Fig. ^, modes can pass through the broad instability regions where many 
of the resonance "tongues" have overlapped. Roughly, the modes lying below A = ^^ 2gj 
will undergo tachyonic resonance, which will dominate the preheating. This matches our 
previous expectation from investigation of the effective mass (|3.6|). 

Our interpretation of excitations of the x field in terms of particles allows us to give 
a physical interpretation of the instability tongues in terms of Feynman diagrams (at least 
in the perturbative regime). The instability regions correspond to those wave- numbers of 
the X field that can be produced by the oscillating inflaton condensate. The only vertex 
available to us is the following: 
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..-•x 



■••••X 

Since we have not quantized the inflaton, we have no inflaton propagators. This allows 
for tree-level diagrams with exactly two outgoing x lines. These diagrams can be classified 
according to the net number of ingoing or outgoing lines of each inflaton field. To associate 
each of these classes of diagrams with a tongue, consider the tongue at (n • m)^ = Am? = 
4A;^a~^. This value represents the energy squared available for the process. To break the 
sign degeneracy, choose n such that n • m > is the energy available for the process. Then 
Hi represents the net number of ingoing (if nj > 0) or outgoing (if nj < 0) inflaton source 
lines. Since two x particles are produced, energy-momentum conservation requires that 
they each have /cphys = k/a = E/2, thus giving E'^ = {n- m)^ = 4A;^a^^. This is precisely 
the location of the tip of the resonance tongue labeled by n. As an example, the first 
diagram corresponding to the the (2,-1) resonance for the case of two fields are shown 
below (where we assume 2mi > m2): 

X 



■ -X 

We will develop the machinery to perform computations in the perturbative regime in 
the next section. The reader should note that we are not claiming the perturbative series 
that can be extracted from our formalism is identical to that familiar from Minkowski 
QFT. Rather, this is intended as a useful picture to understand physically the dynamics 
involved. 

3.7 Four-Legs Interaction: Geometrical Interpretation 

Now, consider the situation for the 4-legs interaction. Defining z = mt, the resonance 
parameters are 

^fc = ^r^ + 2> gi, qi= ^ 3 . (3.17) 

m^a^ ^-^ 4m^a'^ 

A^ now has an explicit dependence on the qis. This results in the restriction 

2^g»<A (3.18) 
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for the 4-legs interaction. Compare this to 2q < A, which would be obtained in the single 



field case. In Fig. 14 where we plot the trajectory through parameter space which is induced 
due to the expansion of the universe. As in the above 3-legs case, we have an identification 
of the various instability tongues with various perturbative decay channels of the inflaton. 




Figure 14: Paths followed by the resonance parameters for several mode functions for the case of a 
4-legs interaction between the inflatons and scalar field. Here we have chosen the parameter q to 
take the value go = 2 at the start of reheating. All of the modes lie above the path foiled by the 
k — mode (orange dotted line). This restriction prevents the modes from passing through the 
broad instability regions. For illustrative purposes, the line A ~ 2q is also shown. This is the line 
the modes must lie above in the single field case. 



Condition ( 3.18| ) is ultimately the reason why nonperturbative decays of the infla- 
ton are absent for the 4-legs interaction. To see this, first suppose that the mass spec- 
trum is such that the condition 2 ^^ Qi 1^ A ensures that the resonances are all distinct.^ 
If we have a large number of g^'s with similar values, then each one satisfies 2qi <^ A 
for ( |3.18| ) to hold. Now consider the value of qi inside the zth first order resonance. We 
have 2(/j <^ A = {mi/m)^, which gives q^ = qi{m/mi)^ <^ 1. Therefore, if the reso- 
nance bands are distinguishable, we have q^^ <C 1 and the decays of the inflatons will be 
perturbative. 

The above argument relies on a notion of distinct resonance bands for each oscillator. 
Let's again assume that all of our squared mass splittings and gj's are equal. From ( 3.1^ ) 
and ( 3.13| ) with A = {rrii/m)^ in the ith resonance band, we see that for none of the 
resonance bands to have overlapped we require 



Nn^r.Sm'^ > 



m,- 



(3.19) 



for all masses mj. Actually, the above estimate is quite conservative. Since qf^ <^ 1 
for each individual oscillator, in order to have a strong resonance we must have many 



*We'll address this point below. 
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different oscillators contributing. Thus, rather than having only pairs of instability bands 
overlapping, we must have a large number of them overlapping simultaneously, which will 
require an even more degenerate mass spectrum than the above estimate. In agreement with 
our investigation of the dynamics of the effective mass ( |3.5| ) , we conclude that provided the 
mass spectrum satisfies ( 3.19| ) the decays via the 4- legs coupling will always be perturbative. 
In this respect we differ from the conclusions of | |53| , |5^ , where the addition of an additional 
oscillating scalar was predicted to increase the efficiency of the resonance. As well, for the 
case of many oscillating fields, the resulting dynamics should be similar to the case where 



the driving term is random [51, 52 1. However, even for the case of many fields, we find here 
that for small q the resonance effects will be too weak to lead to relevant instabilities for 
all but a small subset of modes. This will be confirmed also in section g below. 



To reinforce the point above, in Fig. |15| we show the Floquet exponents along the lines 
A = 2q and A = Aq both for the two oscillator case above and the Mathieu equations of the 
individual oscillators. Although the presence of the second field strengthens the resonance 
along A = 2g, it simultaneously forces the modes above A = Aq actually decreasing the 
strength of the resonance as seen by the mode functions. However, preheating is still 
possible for the case of two oscillators, although it is suppressed relative to the single field 
case Isl 
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Figure 15: Floquet exponents along the curves A = Aq (blue curve) and A = 2q (red curve) for 
both the QP Mathieu (left panel) and periodic Mathieu (right panel) equations. The resonance is 
clearly much stronger along the line A = 2q in both cases. Additionally, the resonance is stronger 
in the QP case along both lines (note the scales on the vertical axes). 



3.8 On the Inclusion of Off-Diagonal Interactions 



In our analysis thus far, we have neglected the possibility of off-diagonal 4-leg interactions 
of the form 

-oflf-diag _ V^ "-'I 
'2 



c 



int 



E^^'^.^' 



(3.20) 



«7^i 



in the prototype action (3T). However, if dj 7^ then no symmetry forbids such terms and, 
in general, they must be included. Let us now argue that their consistent inclusion does 
not add anything qualitatively new to our previous discussion. Note that the interaction 
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( 3.20 ) will yield the following contribution to the effective mass for the x fluctuations: 



AM2ff(t) = ^/i,,</.,(t)0,(t). (3.21) 



Using the solution (2.6) we have 



AM|j(t) = Y^ h,j ^°yy sin [mit + Oi] sin [mjt + % 



«7^i 



where on the last line we have used some elementary trigonometric identities and deflned 



m,^ = rrii i: rrij, 9^- = 0j it 9j + 7r/2. Equation ( 3.22| ) has precisely the same form as the 



t] « JT I] 



contribution from 3-leg interactions of the form ^^ ai(l)iX^. The only difference is the factor 
a~^, which would have been a~^'^ for trilinear couplings. This difference does not play 
an important role for the short time scales relevant for nonperturbative preheating. We 
conclude that off-diagonal 4- leg interactions ( 3.20| ) may lead to nonperturbative preheating 
by tachyonic resonance. 

4. Perturbative Decays of a Single Oscillating Inflaton 

In the last section we studied nonperturbative preheating effects at the end of N-flation. 
In the case of 4-leg interactions we showed that strong nonperturbative effects are absent. 
Hence, the decay of the N-flatons in this case must proceed via perturbative preheating 
decays. On the other hand, in the case of 3-leg interactions we found that strong non- 
perturbative effects are possible via tachyonic resonance. However, perturbative processes 
may still be relevant in this case also, either during the final stages of reheating or when 
the amplitude of the inflaton oscillations is small. Thus, we now turn our attention to 
studying the perturbative decays of the homogeneous oscillating inflatons after N-flation. 
Before we consider the case at hand, A^ ^ 1, it is interesting and instructive to first 
treat the simpler case of single field infiation, A^ = 1, as a warm-up exercise. The material 



in this section is partially review |58, 59 1 and hence the reader who is already familiar 
with the theory of reheating may wish to skip ahead to the next section. Here we explain 
the failure of the infiaton condensate to decay via 4-leg interactions and also show that 
3-leg interactions do allow for a complete decay. Later, we will show that the same results 
persist in the interesting case with N ^ 1. 

As in our multi-field analysis, the homogeneous oscillations of the infiaton are described 

by 

</>(*) = ^5^ sin(mt + 0). (4.1) 

The universe expands as a{t) ~ t'^'^ and the energy density in cp red-shifts like non- 
relativistic dust: p^f, = ^cfP' + '^cfP' ~ a~^. As in ( ^ ) we couple the infiaton to a single 
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scalar field via both 4-leg and 3-leg interactions as 

.2 1 



jr = -^d^^d^^cP-'^^ 



Tidf^xd^X 



5 ,2 o- 
2 ^ 2^ 



2 _ ^ 4 

X ^x ■ 



(4.2) 



The term Ax is required to keep the potential bounded from below, but otherwise will 
play no role in our discussion. We restrict to A > a'^ /{2m?) so that the potential has a 
single minimum (j) = x = with ^ = 0. From (^]^) we can derive the equation of motion 
for the linear inhomogeneous fluctuations of the x field in the background of a classical, 
homogeneous inflaton condensate (p: 



V 



X + 3HX-—X+ [g^<P\t) + acPit)] X = . 



(4.3) 



■^X^) interactions in (4.2) may 



Notice that the combination of 3-leg (</>x^) and 4-leg 
arise in a number of ways. For example, if the inflaton potential V{(j)) has a minimum at 



V about which V^cp) 



^ m2( 



v) /2 then, after the shift 



-|- f , we have the 



familiar form V = m?(j)'^/2. A bi-linear coupling such as g^fp'^x^ then leads to the addition 
of a 3-legs interaction acpx^ (with a = 2g'^v) after such a spontaneous symmetry breaking. 



The scalar potential in (O) may also arise from the super-potential 

9 



W 



m 



+ 



2\/2 2\/2 



:<^X^ 



(4.4) 



which gives A = g'^/2 and a = gm. In our perturbative analysis below it will not be 
necessary to assume that either the 3-leg or 4-leg interaction dominates, we can treat both 
simultaneously. 

4.1 Bogoliubov Calculation 

We now study the perturbative production of x particles from the decay of the oscillating 
fleld ([4.1[). We treat the inflaton field as a classical background and study the linear 
quantum fiuctuations of the inhomogeneous fleld x(^)X). It is convenient to work with 
conformal time r, deflned in terms of the usual cosmic time variable t as adr = dt, and 
introduce a re-scaled "co-moving" field 

/(r,x)=a(r)x(r,x). (4.5) 

We decompose the co- moving field into q- number annihilation/creation operators ak, a^ 
and c- number mode functions /^(r) as 

d^k 



/(r,x) 



a^fk{r)e^'''' + aUt{T)e- 



(27r)3/2 L 
The annihilation/creation operators obey the usual commutation relation 

[ak,ak']=<5(3)(k-k') 
and the mode functions satisfy the oscillator-like equation 

mr)+col{T)MT)=0, 



(4.6) 

(4.7) 
(4.8) 
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where the prime denotes a derivative with respect to cosmic time /' = drf and the time- 
varying frequency is 



u^l{T) = k' + a'[g'c^\r)+aHr)]-^ 



a 
e + ^ sin2 [mt{T) + 6]+ 0^/2(^)^00 gin [mt(r) + 6] , 



a{T) 



(4.9) 



where i(r) = J adr. On the second hne of (O) we have used ( [4.1| ) and the fact that 
a" /a ~ a'^H'^ is neghgible compared to the modes k"^ ~ a'^m? that we expect to be produced. 
(Recall that H -^ m during reheating.) 

We seek a solutions of equation ( [4.8D of the WKB form 



(4.10) 



-i f UJk{r')dT' +ijLUk{T')dT' 

fk{r) = ak{T)^j==^ + /3,(r) 



\/2wfc(r) \/2u;fe(r) 

--.(-)/^e-/-(-')'^^' + ^/3.(r)y^e+^/-(-')^^' , (4.11) 

where the Bogoliubov coefficients are normalized as |afcp — |/3fcp = 1. The ansatz ( [4.10 , 01] ) 



f'kir) 



affords a solution of (4.8) provided akir), (^^{t) obey the coupled equations 



4W = ;^e+2^/-^(-')'^-'A(r), /3Ur) = ^^e'^^ ^ -^'^^'^'^' a,(r) . (4.12) 

The occupation number n^ is defined by the energy of the mode ||/fcP + ^^l\fk\'^ divided 
by the frequency of that mode 






(4.13) 



where the —-^ comes from extracting the zero point energy of the harmonic oscillator and 
on the second line we have assumed the WKB solution ([4.10| ). The occupation number 
depends only on the negative frequency Bogoliubov coefficient, /3/c- 

Demanding that the field be in the adiabatic vacuum in the asymptotic past gives the 
initial condition Ok = 1, /3fc = 0. In the perturbative regime particle occupation numbers 
remain small |/3fc| <^ 1 so we can iterate ([4.12]) to obtain 



A(t) 



-oo 2a;fc(r') 



-2i 



dT'uJkiT") 



(4.14) 



To evaluate the integral ([4.14|) we consider the limit 



A;2 

— >5 



2j,2 



-TT > (Tfi 



(4.15) 



For the modes k ~ am that we expect to be produced, these limits are equivalent to 
assuming that the resonance parameters of the Mathieu equation are small: qs = acpo/m <C 
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1 and §4 = g^(pQ/m'^ <^ 1. Hence, the limits ( 4.15| ) are compatible with the assumption 
that nonperturbative preheating processes are negligible. 

Given ( 4.15| ) we can approximate e~'^J^kdT ^ ^-ijkdr 
for the frequency ( |4.9| ) we have 

-,2 



= e ^^'^. On the other hand, 



u 



^-i[V* + .]5 



Uk 2ujI ^ '"'"'' ^ dr 



i^l''<t>. 



om 



2^1 



29' 



7.3/2 



■ sin [mt{T) + 9] + a 



cos [mt(r) + 9] 



(4.16) 



where we neglected terms with derivatives of the scale factor, a. Collecting everything 
together, we have 



Pkiv) 



where we have defined 



+ 



a(l)Qm 



dr' ( e'^'^4-ic. 



Xr') 



«*:V'4_ieg(^') 



^^/^3/2 ^^iki>U.^(r') 



t 



± 



-2t± 



2mt{T) + 29 
k 



and 



i^tUr) 



-ikip^ 



-2t± 



Ar') 



mt{T) + 9 
k 



(4.17) 



(4.18) 



The first term in ( 4.17 ) is the contribution from the 4-leg interaction (j)(f) — )■ xx while the 
second term is the contribution from the 3-leg interaction (j) — )■ xX- Notice that the integrals 
appearing in ([4.17|) are defined by the interference of two oscillatory terms, e'™*^"^ '~^^ and 



-2ikT' 



Qj, g2«mt(r )+26» g^^^ g-2ifcr ^ -^j-^ere the function t{T) = J adr is non-linear. We will 
give a physical interpretation of these factors later. 

The stationary phase approximation tells us that the integrals appearing in ( 4.17] ) are 
dominated by the contribution near the instant where 



^V'^leg(-) = . 



(4.19) 



For the 4-leg interaction, this is the moment r = T4^fc when 

ma(r4,fc) = k , (4.20) 

while for the 3-leg interaction the dominant contribution to ( 4.17] ) comes from the moment 
T = Ts^k when 

ma{Ts^k) = 2k . (4.21) 

These results are simple to understand on physical grounds: they reflect energy conserva- 
tion for an annihilation of two inflaton particles at rest into a pair of x's, and decay of a 
single inflaton into a pair of x's, respectively. 

Now we can calculate the co- moving occupation number ( 4.13) ): 

|2 (4.22) 

4 ,4.^3 



nk = |/3fe| 



TTCr 



2j,2 



Trg (pn.ni „ r . , , „ r , 



G \a{t) - a-i^k] [a3,fc - ao] 



irg'^acpQ-m^/'^ 



16A;9/2(/74,,/f3,,)i/2 



sm 



V'^-leg(*4,fc) - ^:[-leg(*3,fc) © Ht) - 03,^] 9 [04,^ - a{t)] . 
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The formula ( [4.22| ) is the main result of this subsection. The subscript 4^/; indicates that 
the quantity is to be evaluated at the moment r = T^^k, when k = ma and the 4-leg 
interaction proceeds for a given k. Similarly, the subscript 3^^ indicates that a quantity 
must be evaluated at r = T3 ^ when 2k = ma and the 3-leg interaction proceeds for a given 
k. We denote the scale factor at the beginning of reheating by ag. The step functions 
appearing in ( 4.22| ) enforce the fact that in order for the x's to be produced perturbatively 
in pair annihilation and decays of inflatons by time t, energy conservation requires that 
mao < k < ma{t) and ?n,ao <2k < ma{t) respectively. 



4.2 Boltzmann Equation 

We define the co-moving energy density in the x field as 
d^k 



4 
a p^ 



(2vr)' 



■ojkrik 



9'4 
1287r 



ma{i) 



dk 



m" 



+ 



327r 



H4,kk^ 

ma(t)/2 



+ 



a 



2^2 



327r 



ma{t)/2 



dk- 



1 



mao/ 2 



dk 



m. 



3/2 



mao 



k^/^y^H^^kHs,k 



sm 






4-leg'' 



{tA,k) - i't-leM,k) . (4.23) 



The time rate of change is then given by 



-4 "^ ^ 4 \ 



+ 



g (p^m a (f)Qm 
1287ra6 ^ Mna^ 

V2g^acj)l I H{a) 



327ra9/2 U H{a/2) 



sm 



V'i.lcg(^(«/2))-V'3+_leg(^(«)) 



(4.24) 
XX and 



The first two terms are the familiar result for the Boltzmann equation with (fx, 
(f) — )■ XX processes contributing to the collision integral. The final term in this expression 
arises due to quantum mechanical interference between the two decay channels of the 
inflaton. It is present because the inflaton is a condensate (which we treat as a classical 
source for the quantum x field) and not a collection of particles. 



Notice that the oscillatory "interference" term in (4.24) disappears when we average 
over time and we are left with the usual Boltzmann equation 



_4 u / 4 N 

a ^(aPx) 



-^^^V\.u=Q 2 



+ r 



m 



P<l, + '- <t>^xxP'i> ' 



(4.25) 



where we have used the fact that the energy density in the inflaton oscillations is p^j, = 

The factor of 2 in front of the annihilation cross-section in (4.25|) 



/)2 + m^d?) 



m'^4>o 



arises because there are two inflatons annihilating in each interaction. 

By comparing ( 4.25| ) to ( 4.24 ) we can identify the annihilation cross-section at zero 
relative velocity 



^XX^\v=0 



9 



QAirm? 



and decay width 



(A-i-XX 



327r?7z 



(4.26) 
(4.27) 
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The rate ( 4.27 ) agrees with the standard result obtained in perturbative QFT for the decay 



of inflaton particles via (j) — ?■ XX- The cross-section (4.26) also agrees with the perturbative 
QFT result for inflaton particles at rest, provided that the Feynman amplitude is evaluated 
at V = 0, so that the two incoming inflatons are treated as identical particles in the 
combinatorial counting. (Treating the inflatons as a classical source term accounts for this 
automatically.) 

As usual, conservation of energy-momentum requires a^^-^{a^ p^) = —a~^-^{a'^p^). 
This allows us to write ( 4.25| ) as 



-(a%) = -2 ^'''t"P^^^'^l -=0i.3^,,2 T^, ,„3, 

dt ma^ 



^(a'p^) = -2 ^ "^-^^3 ^--" (a%)^ - r^^^^(aV^) . (4.28) 



This equation has a simple physical interpretation. In the absence of any interactions, it 
shows that the co-moving energy density of inflaton particles remains constant: a^Pfj, ~ 
const. On the other hand, if we include only the decays (j) ~^ XX then equation ( 4.2^ ) 



shows that the co- moving energy density must decrease exponentially: a^p,p ~ e~'"*. In a 
time of order F"^ the inflaton condensate has decayed completely. 

It is interesting to consider equation ( 4.2^ ) in the special case where annihilations 



(j)(j) — )• XX are present but decays (j) — )■ xX a^'e excluded. In this case we can integrate ( 4.2^ ) 
to show that a'^Ptf, — )• const as t — t- oo. In other words, the inflaton does not completely 
decay and some finite co-moving density of non-relativistic (f) "particles" freezes out. The 
problem is that the volume dilution of the non-relativistic particles proceeds faster than 
the annihilation process (l)4> — ^ XX can drain energy from the condensate. Typically, this 
failure to decay will leave a late-time universe which is cold and unsuitable for life, not 
unlike certain parts of Canada. Clearly this is an unacceptable scenario. Our analytical 
conclusion that 4-leg interactions do not allow the inflaton to decay is consistent with the 



results of numerical lattice field theory simulations, see [30| for example 



4.3 Comparison to Inflaton Particle Decays 

In this section we have studied the decay of the coherent, classical oscillations of the inflaton 
field ( |4.1| ) into quanta of the x field. It is instructive to compare our results to what one 
would have obtained for inflaton particles with the couplings ( |4.2| ) to matter. The key 



discrepancy is the appearance of the "interference" term on the last line of ( 4.24 ), which 
averages to zero. As stated above, this term represents quantum interference between the 
decay channels of the inflaton and it arises because we have a coherent inflaton condensate 
(allowing us to treat cp as a classical source). 

To be clear: in our formalism we are describing the production of pairs of x particles 
from a coherent classical inflaton condensate, as opposed to the production of x's from 
annihilations/decays of inflaton particles. These processes are essentially different, even 
though the final decay products are the same. For the decay of classical inflaton oscillations, 
the amplitudes for the 3-leg and 4-leg interactions must be added coherently. This results 
in an interference term when the amplitude is squared - the sinusoidal term on the last 
line of ( PI) . 
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5. Perturbative Inflaton Decays via Four-Leg Interactions 

5.1 Bogoliubov Calculation 

In section ^ we excluded nonperturbative effects for the 4-legs interaction J^idf^lx^- We 
now study perturbative decays of the inflatons via this couphng. Although we have the 
specific example of N-flation in mind, the methods in this and the following section should 
apply to the late time perturbative decay of the inflaton condensate in any multi-field model 
where the condensate oscillates around a minimum in its potential with small amplitude at 



late times. Our approach follows very closely the formalism developed in subsections ^A 



and |4.2| . The reader is encouraged to review these sections for definitions and more details 
of the derivation. The quantum eigenmodes of the xi't,^) field obey the oscillator-like 
equation (|4.8| ) with effective frequency 

^-^ a 

i 

a[T) ^ — ' 

i 

On the last line we have used (|2.6| ) and again neglected a" /a ~ a^H^. The WKB form 
(4.10) still holds with a^ = 1 and /3fc given by the expression ( [4.14 ). 



We wish to evaluate the integrals in (4.14) for the relevant effective frequency (pj). 



To this end, we can consider the short wavelength limit 

i 

To see that this is the relevant assumption, note that modes with k ~ rriia are expected 
to be produced. Therefore the condition (|5.2| ) corresponds to the smallness of the effective 
resonance parameter in the Mathieu equation picture: X^i 9i ^ ^i ^^^ subsection |3.4| . 
Hence, limit ( ^.2| ) is compatible with our previous claim that strong nonperturbative effects 
are absent. 

Using our background inflaton solutions ( p. 61) and keeping only the leading order terms 
in H/rrii, ( [4.14 ) becomes 

/5^(^)-ipE4^^(-)' (5-3) 



where we have defined 



rW, 



atKimi r 



, 2mit{T) + 29i 
i^i KV = -2r ± . (5.5) 

The integral is dominated by the contribution near the instant r = Tj ^ when 

■niia{Ti^k) = k . (5.6) 
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Figure 16: The time-dependence of the co-moving occupation number Uk of x particles with mo- 
mentum k, produced via perturbative 4-legs interactions. We have plotted both the approximate 



analytical estimate (5/7) (the dashed curve) and also the result of an exact numerical computation 
(the solid curve) of the co- moving mode functions fk{T). We take iV = 20 inflaton fields with 
masses ranging from O.GItti to IOtti. For simplicity, we assume equal initial amplitudes 0o,i — 4>a at 
the onset of reheating and also equal couplings gf — "^ ^^° ^^^^^^ic^^a ^^^ r.'^A.'i , 
k'^/m'^ = 1000. 



g . The parameters are g (p^/m = 0.05 and 



As before, (|5.(]| ) has a simple interpretation: it reflects energy conservation for the annihi- 
lation of a pair of (pi at rest into a pair of massless x particles with physical momentum 
k/a. 

(i) 

We can organize the inflatons so that m,j < rrij for i < j. Upon evaluating J^ , the 
expectation value of the co-moving occupation number of the x particles is given by 



nic = |/3fc 



vr 



64/fe3 Z^ 



+ 2E 



j>i 



4 J.4 

9i (f>o,i 



9h'j4>l,i4>l,j 



Q{k,k — to) 



(5.7) 



'i,fc 



Hi^kO-^j^Hjk 



cos{k[ij+{ti^k) - i^^ihk)]) Qihk - to) }@it - ti,k) , 



with '4>f{T) given by (^), ti,fc the cosmic time corresponding to conformal time Ti^k and 
to corresponds to the onset of reheating. We illustrate the accuracy of this approximation 
in Figure |l^. This is an accurate approximation provided |/3fc| <C 1. 
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5.2 Boltzmann Equation 

The energy density is 



rm,a{t) „2 i,2 2 i,2 



-%-T^\T.I d'^^^ (5-8) 






The rate of change of the radiation energy density is 

""'^ dt ^"^^^^ ~ 128^ ^ { ^^^'^0.^"'^ (^-^^ 

As in the single field case, the cross-terms in ( |5.9D arise from QM interference between 
the decay channels. However, unlike section ^, in this case the interference it is due to 
coherence between the same decay channel (i.e. pair annihilations) for different inflaton 
fields. On the other hand, in the single field case the interference was between two different 
decay channels for the same inflaton. 

Assuming a non-degenerate mass spectrum and averaging over a sufficiently long time 
interval for the interference term to vanish, we find 

a-'j^ia^P,) ^2j2^^M^^f^pl , (5.10) 

i 

where we have used the fact that the energy density of the i-th inflaton is given by 



2 o ,ox - "*: 



2j,2 



^{(pi +mf(l)f) = ^'^3'° . By comparison of ( |5.1C1| ) to (5^) we can extract the pair 



annihilation cross section for two 4>i^s into a pair of x particles at zero relative velocity is 

K.^.^xxvl=o = ^ ■ (5-11) 

This is consistent with the perturbative QFT result one would obtain for pair annihilations 
of inflaton particles, provided one accounts for the subtlety concerning the limit w — )• 
discussed in section ^. 

There is an interesting subtlety associated with the time-averaging that leads to ( 5.10| ). 



For mass spectra which are very nearly degenerate {i.e. m^^ — m^jj^ <C "lavg) the time 
over which we must average may be much longer than m^ and at any given time the 
oscillating terms in (]^M) may contribute more than the non-oscillating terms. This is not 
surprising, since over time intervals much less than the inverse mass splittings between the 
various inflatons, we may expect the many individual inflatons to behave as a single inflaton 
with oscillation amplitude 4> ~ ^i't'i- In this case, our effective resonance parameter will 
be ^^ qi where the sum is over those inflatons whose masses are nearly degenerate, at least 
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at early times before the inflatons have a chance to de-phase. If almost all of the inflatons 
are degenerate in mass, it is possible for this effective resonance parameter to become large 
and the resonance to be efficient. As a simple example, if all the inflatons have equal mass 
and the same initial conditions, then all of the cosines will be 1 and the "oscillating" terms 
(which there are N"^ — N of) will contribute more than the "non-oscillating" terms (which 
there are A^ of). Note that when we say degenerate, we mean that the mass splittings 
between all pairs of inflatons are small, not just the mass splittings between rm and rriij^i. 
Assuming that the quantum mechanical interference terms are negligible, we can pro- 
ceed to write down the effective Boltzmann equation for the inflaton oscillations. Conser- 
vation of energy-momentum gives us 

«^'|PE/'^»)+«"'1(«'/'x) = 0- (5-12) 

From ( 5.12| ) and ( |5.1C1| ) we have the Boltzmann equation 

37 2^ "Z'* J =-2 2^ -—3 {a P<^J ■ (5.13) 



By direct integration one may verify that 



a^ 



^P0. -> const (5-14) 



as t — )• oo. Hence, the decay of the N-flatons does not complete: some finite co- moving 
energy density always freezes in. This is quite analogous to the result for single field 
inflation with 4-legs interactions, see section ^. We conclude that the model ( |3.lD with 
(Tj = is not, in general, a viable scenario. 

6. Perturbative Inflaton Decays via Three-Leg Interactions 

The 4-legs interactions discussed in the last section do not permit a complete decay of the 
inflaton. Therefore we must include other couplings between (pi and matter. A natural 
candidate is the tri-linear (3-legs) interaction (piX^. In this section we develop the theory 
of perturbative reheating after inflation for the model ( |3.1| ) with gf = but ai ^ 0. For 
the nonperturbative decays, see section ^ 

6.1 Bogoliubov Calculation 

Let us now discuss the perturbative decays of the inflatons in the presence of 3-leg interac- 
tions J2i '^i4'ix'^- Such perturbative decays are relevant during the final stages of reheating, 
or in the case where the inflaton oscillations at the end of inflaton are small. The quan- 
tum eigenmodes of the x(^) x) field obey the oscillator- like equation ( [1.8| ) with effective 
frequency 

a" 



a 



^^ + 1/9/ N V! <7j(/>o,j sin [mit{T) + 6i\ , (6.1) 
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where on the last hne we have used (|2.6| ) and again neglected the a!' ja ~ a^-ff^ term. 

Following our previous formalism, we wish to evaluate the integrals in (^.14) for the 



relevant effective frequency ( |6.lD . The dominant contribution to the integral for the i-th 
inflaton comes from the moment when 2/c = amj, corresponding to energy conservation for 
the process 0i — )■ XX- We find the following result for the occupation number: 



nk = |/3fcP = -^5] 



cr; 



2x2 



0(aj,fc - ao) (6.2) 

0(a-ai,fe), 



+ 2 V] -^ J" ^ '^ cos(V'i(ti,fc) - i)j{tj^k))e{aj^k - ao 



where ^i(t) = — 2A;r(t) + rriit + Oi, rriiai^k = fniO'{ti,k) = 2fc, ffj^fc = H{ti^k) and m^ < ttij 
for i < j. 

The accuracy of this analytic calculation is illustrated in Figure ^. The oscillations 
appearing on the graph arise from boundary terms that have been ignored in the stationary 
phase approximation (i.e. higher order terms in the complete asymptotic expansion of the 
phase integral ( [1.14D ). Physically, the variation of the background value of the inflaton 
field leads to continual creation and destruction of x particles. However, except in the 
vicinity of the stationary points 2k = rriia, the destruction of particles tends to periodically 
balance the creation of particles, leading to the oscillatory behaviour superimposed onto 
our analytic solution. 

6.2 Boltzmann Equation 

We find the average rate of radiation energy production 






^: («Vx) = ^Y1 \Z ' =^Y1 ^'t'^^xxP'P. . (6-3) 

where 



I I 



<y\ 



Here we recognize V^^^^^ as the decay width for c/ij — )• XX due to the coupling ^Oi^iX^ . 
This coincides with the result one would compute in perturbation theory for the decay of 
inflaton particles. 

Conservation of energy gives a~'^dt{a'^p^) = —0'~^J2i'^t{o-^P(f>i), hence we can write 
( |6.3| ) as a series of N decoupled equations 

^(aV<^J = -r</,,^xx(a^/'<AJ ■ (6-5) 

By direct integration we see that the co-moving energy density of each inflaton drops 
exponentially a^p,fn ~ e~'"'* and hence the decay process does complete. This is similar 
to the single field case with 3-legs interactions, but unlike the case with 4-legs interaction 
(both single and multi-field versions). 
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Figure 17: The time-dependence of the co- moving occupation number Uk of x particles with mo- 
mentum k, produced via perturbative tri-hnear interactions. We have plotted both the approximate 
analytical estimate (|6.2| ) (the dashed curve) and also the result of an exact numerical computation 
(the solid curve) of the co- moving mode functions fk{T). We take iV = 20 inflaton fields with 
masses ranging from O.OItti to IOtti. For simplicity, we assume equal initial amplitudes 0o,j — 4>o 
at the onset of reheating and also equal couplings ai = a. The parameters are aipo/m^ = 0.05 and 
k^/iTi^ = 1000. Compared with the 4-legs interaction in Fig. HO, the particle occupation numbers 

(fc/m)3/2 



are greater by a factor of a.^j. 



10° as expected from comparison of (5.7) and (3.2) 



7. Reheating into the Standard Model 

7.1 Coupling to Gauge Fields 

The prototype Lagrangian (^]^) allows us to study the decay of the inflaton condensate 
into scalar field fluctuations from a phenomenological perspective. However, in order to 
truly make contact with the usual picture of the hot big bang, we would like to understand 
in detail the production of standard model degrees of freedom at the end of inflation. A 
novel feature of string theory inflation model-building is that it is possible, at least in 
principle, to determine such couplings from a bottom-up perspective. As a first step to 
understanding reheating of the SM after N-flation, we consider a U{1) gauge field A^ with 
field strength F^i, = d^A^, — d^A^. Since (pi are axions, we expect couplings to F A F. 
Hence, we consider the action 



N 



^ = E 



i=l 



-i^di,(t>id^(t>t - ^^ 



If, 



liv 



p^iu 



E 



4M,; 



-F'^^F, 



^IV 1 



(7.1) 



where F^'^ = 2^^'^^'^Ppa- is the dual field strength and Mj is some symmetry breaking scale. 
Typically Mi ~ /j, however, we do not require this for our analysis. 

It is worth considering how ( [7.lD arises in the IIB string theory model discussed above. 
Since the axions (pi are associated with 4-cycles, it is natural to consider localizing the SM 
on a D7-brane wrapping the j-th cycle. (See also |11, ^.) In this case one expects a 
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coupling of the form 

Ant = -X^jFF 

to brane-bound gauge fields. In general, the mass basis 



(7.2) 



bi will not be aligned with the 
basis (^j associated with the geometrical 4-cycles. Hence, a rotation in field space ( ^.^ ) 
will be necessary to put the axion action (2.7) in the canonical form ( |7. 1| ) . For a general 
rotation all the mass basis axions 0j will appear in the sum (^j^) and, accordingly, all of 
the N-flatons are expected to couple to F A F. This argument leads directly to an action of 



the form (7.1) where the couplings M, are related to the coefficients oj- in the field-space 



.W 



1 



rotation ( p.8[ ). These are, obviously, very model dependent. Normalizing ^ • 

we expect af ~ N^'^/'^ so that Mi ~ N^/'^/X. 

Note that the argument above works equally well for gauge fields living on any D7- 
brane wrapping any 4-cycle of the Calabi-Yau compactification. In general, there might 
be many such branes so that the compactification contains the usual SM in addition to 
a large number of hidden sectors. Those axions which are still dynamical at the end of 
inflation will reheat into all sectors more-or-less equally (of course, in a specific model this 
depends on the details of the field rotation and the placement of branes in the Calabi-Yau 
manifold). Relics in the hidden sectors might be dark matter candidates, or they might 
pose cosmological problems, depending on model parameters. See I^Tj for more details. 

is 



The equation of motion for the gauge field derived from (7.1 



i)F^'•' = 0. 



(7.3) 



As usual, we study the linear quantum fluctuations of the matter field, treating the inflaton 
oscillations (2^) as a classical background. We work in a transverse gauge, with ^"^ = 
and diA^ = 0, and decompose into circularly polarized mode functions 

A\t, x) = I j0j^ Yl [e\(k)a^e^'^-yl^(t) + h.c] , (7.4) 

where the circular polarization vectors satisfy k x ej- = ^ikepm, o^ ^^^ q-number anni- 
hilation/creation operators and "h.c." denotes the Hermitian conjugate of the preceding 
term. Using this expansion, (^^) gives 



A^ + HA^ + 



k^ k^ 
~^^~a^ 



m 



M, 



A^ 



0. 



(7.5) 



7.2 Nonperturbative Preheating 

To investigate the possibility of nonperturbative photon production, we would like to put 
(|7.5|) in the form of the QP Mathieu equation. To this end we introduce the co-moving 



modes A^ = a^''^A^ and define a new time variable 2z = mt where m is a typical inflaton 
mass. Treating the expansion of the universe adiabatically one arrives at the QP Mathieu 
form ( 3.1C| ) with parameters 

4A;^ I 2kmj(bnj 



m 



2M,-a5/2 



(7.6) 
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A=10q^ - 



A=0.1q^ — 
A=q^ — 

Figure 18: The permitted regions for gauge field modes in tlie (A, q) plane. The field modes must 
lie above the curves A = C^'^{2q)'^, where the parameter C = J2i 'fni(t>o,i/mMi. 



We can now employ our results from subsection 3.4. Notice that, unlike the scalar field 
cases studied in section ^, both our Aj. and qi parameters now depend on k. The resonance 
parameters for the mode functions are no longer initially distributed along a vertical line, 
but instead along the parabola 



Al 



m Mi 



W) = E 



rriiCpQA 




mi (j)Q^iJ \'^-^ mMi 

The expansion of the universe causes these parameters to "flow" along the lines 



(7.7) 



Ak 



Al 



4/5 



Ak 






4/5 



(7.8) 



Hence, all of the field modes will trace out curves in the {A, q) plane that lie above 
(X^i 9«)^ — C!^^iJ2i li)"^- This is illustrated in Fig. 18 for several choices 



A 



y.i 



% m Mi 



of the prefactor C. In Fig. [1^ we show the Floquet exponents along several of these initial 
curves, which provides a measure of the initial instability of the gauge fields. 

For Qi values such that the individual bands are distinct, the condition qf ^ 1 while 
the mode is in the ith first order resonance band at A = {mi/m)'^ translates into 

^)' « ' ■ <"> 

90,1 J 

This is the straightforward generalization of the condition found in |5^. Meanwhile, in 
order to for the parabola to pass through the region where many of the first order bands 
have overlapped and we again expect the resonance to be strong, we instead require 



Ndrn^ 



E 



mz 



NoscSm'^Mi 



> 1 



(7.10) 
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Figure 19: Floquet exponents along the curves A = C ^{2q)^ for several choices of C and the same 
parameters as in Figure |l^. 



where we have again assumed that all of the squared mass splittings 6m'^ and qi = g's are 
equal. 

It is worth considering whether the condition for an efficient resonance is satisfied 
in the IIB string theory construction discussed in section g. As mentioned above, we 
expect a coupling (j)jFF/M to brane-bound fields living on the j-th cycle. In string theory, 
we typically have M below the Planck scale, but not significantly so. Thus, we assume 
M ~ Mp. Upon rotating to the interaction basis we have a coupling to all of the N-flatons 

N"^/"^ (as discussed above) which leads to 



(j)i. For a generic rotation (2 
the estimate M, 



we have ai 



NMp for the couplings in ( |7.1| ). Combining this with the estimate 
osc gives (pi^o/Mi ~ \/^JWNo^c ~ VlQ/N . In this case the condition ([Tsl) 



for a single one of the fields to provide an efficient resonance is not satisfied. Meanwhile, in 
order for us to be in the regime ( |7.10| ) where many resonances have overlapped we require 



"^max^^ifAci 



\ ^ ' ' ''max 



-'*osc'^typical'^max 

N5m? 



m. 



typical 



>1. 



(7.11) 



Taken together, these constraints imply that, as long as the mass spectrum is not too 
degerate, nonperturbative effects will be weak and the decay of the inflaton will be domi- 
nated by perturbative processes (pi — )• 77. Obviously, this conclusion will depend consid- 
erably on the model-building detail, especially the field rotation ( |2.8D and the presence of 
hidden sectors in the Calabi-Yau compactification volume. 



7.3 Perturbative Reheating 

With the expectation for perturbative decays described above, we now investigate this case. 
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Transforming ( |7.5D to conformal time gives 

d2 



We see that our modes functions A^^ are already co-moving with an effective frequency 
given by 



{, ,±^2 _ i,2 _ ^ V^ 4>0,imi 



cos{mit + 9i) — - — sin(mjt + 6i 



(7.13) 



The second term in the square brackets is small compared to the first term during reheating 
and can be ignored. Following our procedure from before (replacing fk with Aj^ , and 
assuming ^^ ^' "Q- <C 1), we find for the co-moving number density (ignoring the cross 
terms) 



,± 



/ \ -^/^ 5/2 

vr ^^ mfcl)l- vr^/12 1 /^ 2.2 1 ^^ ^0,iW 



ITTm = ^^Tm Z^^^^'^o,, 2^ — TT^ , (7.14) 



J , , J 



where we have assumed that the energy density of the universe is dominated by the inflatons 
{i.e.H"^ = Q^3j\^i ^i ^00 i) in the last step. Accounting for both polarizations and time 
averaging, the gauge field gains energy at a rate 

i ^ i 

The decay rate for (pi is given by 

3 

This agrees with the more familiar result for the perturbative decay of inflaton particles. 
As before, we can use conservation of energy to re- write ( 7.15| ) as 



so that a^Pcf)^ ~ e~'"** — )• as t — )• oo, indicating that the N-flatons are permitted to decay 
completely into photons. 

Let us now estimate the reheat temperature provided by the above decays. When 
decaying into a single field, we have the standard estimate for the reheat temperature 

Tr^{-^\ ^/MpT (7.18) 

where 

r = yri^ (7.19) 
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and g is the number of relativistic degrees of freedom (including A) which interact with 
our "photons". Let's assume that the inflatons ah have similar masses and oscillation 
amplitudes. Then P(j,jp(f, ~ N~J. and we have 



r«H^V",/^y^^^~J^~i5!^. (T.20) 



gny \64tt^NoscM^ V NMp ,/n 



where we used rrii ~ lO^^Mp for Nflation. Notice that the suppression factor is the total 
number of inflatons (which determines the rotation in the field space), not the number of 
inflatons contributing to the reheating. 

8. Applications to More General Models 
8.1 Generalizations 



Clearly, much of our analysis in this paper is not specific to N-flation |18] or its IIB string 



theory realization [19|. The general form of the potential (2.2) will hold near the minimum 
in a large class of models. Thus, the decay of the coherent oscillations ( |2.6D might be 
relevant for a variety of multi-field constructions. In this subsection we speculate on some 
possible applications of our results outside of the context of N-flation.^ We leave a detailed 
discussion to future works. 

In particular, the decay of the inflaton oscillations ( |2.6D via couplings such as ( |3.1| ) 
and ([7.1D might be important for models of multi-field inflaton on the string theory land- 



scape [£l|, ^, ^, ^, ^. The inflationary trajectory across the cosmic landscape typi- 



cally involves a large number of twists, turns and bifurcations, leading to various features 



p3| , plf and nongaussianities |25| in the Cosmic Microwave Background (CMB) fluctua- 
tions. Inflation might end when this trajectory encounters a steep minimum, leading to 
a post-inflationary evolution where are large number of scalars (pi oscillate according to 
( |2.6| ). Depending on the speciflc corner of the landscape, the mass spectrum rrii and initial 
conditions at the onset of reheating might be more-or-less randomly distributed. 



Finally, our analysis might also have some relevance for the model of [42| if the final 



stage of infiation is not entirely along the infiaton direction associated with our visible 
standard model. In that case our results could have relevance for the nonperturbative 



production of hidden sector baryonic dark matter. See also ||61|] for a related discussion. 



8.2 Modular Inflaton Couplings to Gauge Fields 

A more specific stringy scenario where our analysis might prove relevant is racetrack infiaton 
1 43, ^]. In the "improved" construction of P^ j infiation involves the complex motion of 
several moduli and axions in the IIB KKLT vacua discussed above. If there is a prolonged 



^Note, however, that the following discussion will certainly not apply in all models. For example, in 
the presence of steep tachyonic directions the energy of the inflaton will be mostly drained by spinodal 
decomposition before the oscillations (p.q) can set in p] . Also, it may happen that the potential has a very 



steep minimum so that the form (2.2) is not valid for the entire duration of a single oscillation, as in the 
roulette inflation model jll|. 
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phase of oscillations about the minimum^, then our results will apply. We expect couplings 
of the form ( |7.1| ) for the inflatons associated with axions arising from wrapping the 4- 
form field on various 4-cycles of the Calabi-Yau compactification manifold. However, there 
will also be couplings between brane-bound gauge fields and the moduli which control the 
volume of these 4-cycles. Such couplings typically take the form 0j-F^ for gauge fields living 
on the j-th 4-cycle [0]. On rotating to the mass basis ( p.8[ ) we will have couplings between 
many inflatons and the field strength squared. Thus, in addition to the axion coupling 
( |7.1| ), there is a phenomenological motivation to consider "modular" couplings of the form 



N 



^ = E 



i=l 



\d^(p.,d^cl). 



ml 
2 



1 



F,„.F'"' 



flU 



^4M, ^"^ 



Here we assume ^^ (l)i{t)/Mi ^ 1 in order to justify the neglect of higher order terms in the 
infiaton potential and of dimension-6 and higher interaction terms in (|8.1| ). Interactions of 
the type considered in (SJ.) might play a role in a model of assisted inflation from string 
theory using the dynamics of moduli and, perhaps, also in the model of [38|. Working in 
transverse gauge, we find that our mode functions obey 



4" + 



^+E 
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f^2 

4 + ^4 



0. 



Here we work to leading order in ^^ (pi(t)/Mi <^ 1. To put this into the form of a time- 
dependent frequency oscillator, we define a new field A = ^/oH'A. The equation for our 
mode function now becomes 



4 + 



fc^ 
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2'^ M.; 
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iff' 
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A^ 
^k 
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(8.3) 



Again, we drop terms which are higher order in X^j '?^j(^)/-^j ^ 1- The fact that the 
modular coupling in ( ^.l] ) must be treated as a small perturbation to the gauge field 
kinetic term rules out strong nonperturbative preheating effects during the stages when 
the ^j "T-f 0f form of the potential is valid. Hence, we study instead the perturbative 
decays of the coherent infiaton oscillations into gauge field quanta. Performing our standard 
procedure as in section ^ and accounting for both photon polarizations we find the effective 
Boltzmann equation 



with decay rate 



d / 3 N 



-imod 



^T\^'p^.) 



rat 



U'kM} ■ 



5.4) 



(8.5) 



This again matches the standard result for decays of infiaton particles. The interaction (| 
permits the complete decay of the infiaton with a reheat temperature Tr ~ A^~^'^10^'^GeV 



®Such a phase does exist in the case of the original racetrack model [U3| where nonperturbative preheating 
effects are absent in the scalar sector [MSl. 
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9. Conclusions 

In this paper we have studied (p)reheating after inflation in models with a large number 
of inflaton fields. Such models are quite natural from the string theory perspective and 
may predict significant gravitational waves in a regime where the individual inflaton dis- 
placements are sub-Planckian, via the assisted inflation mechanism. We have focused on 
N-flation |18| in particular, however, we expect that our results may apply more gener- 



ally to a variety of multi-field inflation models. In particular, our finding may have some 



relevance for multi- field inflation on the cosmic landscape [53, 24, B5[. 

We have proceeded phenomenologically, considering a variety of possible couplings be- 
tween the inflaton fields and matter. We studied the decay of the N-flatons into scalar field 
matter via both 4-leg and 3-leg interactions. In the former case, we found that reheating 
does not complete (similarly to single field inflation). In the latter case, we found that 
reheating may complete in the perturbative regime, and also that strong nonperturbative 
tachyonic resonance is possible for certain parameters. 

In the case of 3-leg interactions, the fluctuations of the preheat field display exponential 
growth in certain regions of phase space. The structure of the stability/instability bands 
is quite rich and is related to the mathematical properties of the quasi-periodic Mathieu 
equation. We have shown that the presence of multiple frequencies in the effective mass of 
the preheat field leads to interference effects and the dissolution of the stable regions. 

We also studied the decay of the N-flatons to gauge fields. In the original proposal 
p^ the N-flatons are axions, therefore we considered a coupling of the form ^,- (piFF. 
We couple all of the N-flatons to a single U{1) gauge field since the mass basis (j)i and 
interaction basis (pi are not, in general, aligned with one-another. We have discussed how 
this interaction arises in the context of type IIB string theory vacua. We found that 
(p)reheating in this case proceeds very similarly to the 3-legs scalar interaction. Reheating 
can complete perturbatively and strong nonperturbative preheating effects are possible. 

For completeness, we have also considered a coupling ^^ (piP'^ to gauge fields. Such 
terms may arise in string theory models where assisted inflation is realized from the dy- 
namics of moduli flelds. 

During the course of our analysis we have developed an analytical theory of the per- 
turbative decays of coherent, classical inflaton oscillations. This formalism applies much 
more generally than the N-flation models considered here. We have used this method to 
compare the difference between the decay/annihilation of a classical, homogeneous infla- 
ton condensate with the analogous process for inflaton particles. We have seen that these 
two processes are essentially different. In the case of a classical inflaton, the various decay 
channels of (p contribute coherently to the amplitude. Upon squaring this, we see that these 
decay channels can interfere quantum mechanically leading to oscillatory contributions to 
the effective Boltzmann equation. These oscillations disappear upon averaging over a suffi- 
ciently long time interval and we recover the usual results for decay/annihilation of inflaton 
particles at rest. 

There are a number of directions for future studies. It would be interesting to consider a 
rigorous, stabilized stringy embedding of N-flation and identify the location of the standard 
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model and the various decay channels of the inflaton |41|. Using the results of this paper, 
one could assess the viability of such a construction and study the production of exotic 
relics, gravitational waves [U], nongaussianities [15, ^], etc. It would also be interesting 



to study post-inflationary dynamics in racetrack inflation [43, 44 1, or in a more general 
string landscape setting [23, p4[ . 

There does not yet exist a complete theory of preheating after multi-field inflation. The 
result of this paper should be considered as a step in the direction of a more comprehensive 
study. As in the single field case, we expect that a rich spectrum of perturbative and 
nonperturbative QFT effects are possible [Q[, in addition to distinctively stringy processes 
ill, Hi- 
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